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^ ■ ABSTRACT 
CN ■ 

, Context. Present X-ray missions are regularly discovering new X/y-ray sources. The identification of the counterparts of these high- 

' energy sources at other wavelengths is important to determine their nature. In particular, optical observations are an essential tool in 

^ ' the study of X-ray binary populations in our Galaxy. 

' Aims. The main goal of this work is to determine the properties of the optical counterpart to the INTEGRAL source IGR J06074+2205, 

Xf^ , and study its long-term optical variability. Although its nature as a high-mass X-ray binary has been suggested, little is known about 

■ its physical parameters. 

Methods. We have been monitoring IGR J06074+2205 since 2006 in the optical band. We present optical photometric BVRI and 

I I , spectroscopic observations covering the wavelength band 4000-7000 A. The blue spectra allow us to determine the spectral type and 

■^C , luminosity class of the optical companion; the red spectra, together with the photometric magnitudes, were used to derive the colour 

a excess E{B - V) and estimate the distance. 

' Results. We have carried out the first detailed optical study of the massive component in the high-mass X-ray binary IGR 

I J06074+2205. We find that the optical counterpart to IGR J06074-H2205 is a V = 12.3 BO.SVe star located at a distance of -4.5 

^ kpc. The monitoring of the Ha line reveals V/R variability and an overall decline of its equivalent width. The Ho- line has been seen 

I ' to revert from an emission to an absorption profile. We attribute this variability to global changes in the structure of the Be star's 

O ' circumstellar disc which eventually led to the complete loss of the disc. The density perturbation that gives rise to the V/R variability 
vanishes when the disc becomes too small. 
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C/3 . 

^ Key words, stars: individual: IGR J06074+2205, - X-rays: binaries - stars: neutron - stars: binaries close -stars: emission line. Be 

J>. 1 . Introduction tected only one source in the INTEGRAL error circle at a posi- 

l/-) tion of RA = 06h07m26s.62, Dec = 22d05m47s.6 (J2000). The 

^ IGRJ06074-f2205 was discovered by INTEGRAL/JEM-X during 2 x 10-'^ erg cm-^ s' in the energy range 0.3-10 keV. 

public observations of JheCrab_region that took place on 15 and ^he energy spectrum could be fitted with an absorbed power-law 

^ 16 February 2003 ( jChevenez et al. ||2004). The source was de- ^^^^ Nh^(6 + 2)x 10^2 cm ^ and a photon index of 1.3 + 0.8 

• , tected with a flux of ~ 7 mCrab (+2 mCrab) in the energy range ^gQ% confidence errors). The absorbed flux value was nearly 60 

Q , 3-10 keV and 15 mCrab in the range 10-20 keV range on 15 ^jj^es lower than the value of its discovery, which is consistent 

Q February 2003. A day later the flux had decreased to less than 5 ^j^^ jj^e behavior of Be X-ray binaries. 

mCrab. From low-resolution spectroscopic observations obtained us- 

— . Optical spectroscopic observations performed with the jng the G. D . Cassini 1.5m telescope in Loiano observatory, 

MDM 2.4 m telescope of the brightest stars in the field ai'ound Ma setti et al.1 (^6) suggested a B8III optical counterpart to 

V the INTEGRAL position on 27 December 2005 found a Be star iqr J06074-H2205. However, their spectra did not have the 

^ within r of the INTEGRAL position, exhibiting Ha in emis sion spectral resolution to allow a precise classification. In fact, no 

. 5^ with an equivalent width of -6.6 A (Halpern & Tva^ill2005h. Its Galactic Be/X-ray binary contains a Be star with spectral type 

optical and infrared brightness obtained from the USNO A2.0 later than B3, which casts some doubt on the accuracy of this 

and 2MASS catalogs is B = 13.3, R = 12.1, J = 10.49, classification. 

H - 10.19, K - 9.96 mag. Nevertheless, the available data seem to indicate that IGR 
On 2 December 2006 a 5-ks Chandra observation improved J06074-H2205 is a Be/X-ray binary. Be/X-ray binaries ai-e a class 
the accuracy of the X-ray position and allowed the confir- of high-mass X-ray binaries that consist of a Be star and a neu- 
mation of the Be star suggested by Habern & Tyagi (2005) tron star (Ziolkowski 2002). The mass donor in these systems is 
as the correct optical counterpart (Toms ick etah 2006). The a relatively massive (>10 Mq) and fast-rotating (-80% of break- 
Chandra observation also served to discard flie radio source up velocity) star, whose equator is surrounded by a disc formed 
NVSS J060718-H220452, which lies 80" from the reported po- from photospheric plasma ejected by the star. Ha in emission 
sition of the X-ray source (' Poolevlllool iPandev et al.l l2006h . as is typically the dominant feature in the spectra of such stars. In 
a possible radio counterpart to IGR J06074-H2205. Chandra de- fact, the strength of the Balmer lines in general and of Ha in 
particular (whether it has ever been in emission) together with a 
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Table 1. Photometric measurements of the optical counterpart to IGR J06074+2205. 



Date 


JD (2,400,000+) 




B 


V 


R 


/ 




25 October 2007 

26 October 2007 


54399.62 
54400.29 


12.^ 
U.i 


!8±0.02 
! 1+0.03 


12.28+0.03 
12.26+0.03 


11.90+0.02 
11.89+0.04 


11.47+ 
11.43+ 


:0.02 
:0.05 



class of objects. Be stars are also obse rved as single objects, i.e . 
not forming part of a binary system (iPorter & Riviniusll2003l) . 
The equatorial discs are believed to be quasi-Keplerian and sup- 
ported by viscosity (IOkazaldll200lh . The shape and strength of 
the spectral emission lines are useful indicators of the state of 
the disc. Global disc variations include the transition from a Be 
phase, i.e., when the disc is present, to a normal B star phase, i.e., 
when the disc is absent and cyclic V/R changes, i.e., long-term 
cyclic changes in the ratio of the blue to red peaks of a split pro- 
file that are attr ibuted to the pr ecession of a density perturbation 
inside the disc (IOkazakilll991h . 

In this work we present the first detailed study of the optical 
counterpart to the X-ray source IGR J06074-I-2205. Our photo- 
metric observations are the first dedicated observations of the 
source. So far, the reported magnitudes come from various cat- 
alogs. In most cases it is not possible to know the exact date 
of those observations, hence it is not possible to perform vari- 
ability studies. In addition, our monitoring of the Ha line allows 
us to investigate the long-term variability of the circumstellar 
disc. We confirm that IGR J06074+2205 is a Be/X-ray binary. 
How ever, its op tical companion is more massive than suggested 
bv lMasetti et all (f2006D . The main body of the paper corresponds 
to Sect. 3 and 4. In Sect. 3 we present the results of our spectro- 
scopic observations, with emphasis on the variability of the Ha 
line, while in Sect. 4 we discuss the implications of the results. 
Sect. 2 describes the observations and in Sect. 5 we present our 
conclusions. 



2. Observations 

Optical spectroscopic and photometric observations of the op- 
tical counterpart to the INTEGRAL source IGR J06074-H2205 
were obtained from the 1.3m telescope of the Skinakas ob- 
servatory in Crete (Greece) and from the Fred Lawrence 
Whipple Observatory at Mt. Hopkins (Arizona). In addition, 
IGR J06074H-2205 was observed in service time with the Nordic 
Optical telescope (NOT) and the EEV42-40, 2Kx2K chip on 
the night of 25 December 2007 (Grism#16) and on the night 
28 March 2010 (Grism#8). The 1.3 m telescope of the Skinakas 
Observatory (SKI) was equipped with a 2000x800 ISA SITe 
CCD and a 1302 1 mm ' grating, giving a nominal dispersion 
of ~1 A/pixel. We also observed IGR J06074+2205 in queue 
mode with the 1.5-m telescope (FLW) at Mt. Hopkins (Arizona), 
and the FAST-II spectrograph plus FAST3 CCD, a backside- 
illuminated 2688x512 UA STA520A chip with 15yum pixels. 
The observation on 12 September 2007 was obtained with the 
12001/mm grating, while the rest with the 6001/mm grating. 
Spectra of comparison lamps were taken before each exposure 
in order to account for small variations of the wavelength cali- 
bration during the night. To ensure an homogeneous processing 
of the spectra, all of them were normalized with respect to the 
local continuum, which was rectified to unity by employing a 
spline fit. 



The photometric observations were made from the 1.3-m 
telescope of the Skinakas Observatory. IGR J06074H-2205 was 
observed through the Johnson B, V, R and / filters. For the 
photometric observations the telescope was equipped with a 
2048x2048 ANDOR CCD with a 13 5 jum pixel size. Standard 
stars from the Landolt list (iLandoItl l2009l) were used for the 
transformation equations. Reduction of the data was carried out 
in the standard way using the IRAF tools for aperture photome- 
try. 

The photometric magnitudes are given in Table [T] while the 
log of the spectroscopic observations is shown in Table |2l 

3. Results 

3. 1 . The Ha line: evolution of spectral parameters 

Our monitoring of IGR J06074-H2205 reveals that the Ua line 
is highly variable, both in strength and shape. The line always 
shows a double-peaked profiles but the relative intensity of the 
blue (V) over the red (R) peaks varies. 

Table |2] gives the log of the spectroscopic observations 
and some important parameters that resulted from fitting two 
Gaussians to the Ha line profile. Column 6 gives the equivalent 
width of the Ha line (EW(Ha)). The main source of uncertainty 
in the equivalent width stems from the always difficult definition 
of the continuum. The EW(Ha) given in Table |2] correspond to 
the average of twelve measurements, each one from a different 
definition of the continuum and the error is the scatter (standard 
deviation) present in those twelve measurements. 

Column 7 gives the peak separation. Apeak, of the blue and 
red peaks. The peak separation is simply the difference between 
the central wavelength of the red minus the blue peak in velocity 
units (AA/A x c), where c is the speed of light. The errors are 
estimated by propagating the uncertainty in the determination of 
the central wavelength of the gaussian profile used to fit the Ha 
profile. The uncertainty of the best-fit peak wavelengths are at 
68% confidence intervals. The peak separation is related to the 
size of the Be star's disc since 

^disc _ pvsin 

Rt \ Apeak / 

where v sin / is the projected rotational velocity of the B star (v 
is the equatorial rotational velocity and / the inclination of the 
equatorial plane with respect to the observer). 

Column 8 shows the ratio between the core intensity of the 
blue and red humps. The V/R ratio is computed as the logarithm 
of the ratio of the relative fluxes at the blue and red emission 
peak maxima (without prior correction for the underlying contin- 
uum level). Thus negative values indicate a red-dominated peak, 
that is, V < R. 

Figure[T]displays the evolution of the line profiles. V/R vari- 
ability is clearly seen, indicating significant changes in the struc- 
ture of the equatorial disc on timescales of months. In addition. 
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Fig. 1. Evolution of the Ha line profile. V/R variability as well as a long-term decrease of the hne strength is seen. 
Table 2. Log of the spectroscopic observations. 
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Date 



JD (2,400,000+) Telescope 



Wavelength 
coverage (A) 



Resolving 
power' 



EWCHci-) 
(A) 



^peak 

(kms^') 



iog{V/R) 



22-10-2006 


54031.58 


SKI 


5845-7200 


2200 


-11.5+0.4 






24-10-2006 


54033.61 


SKI 


5845-7200 


2200 


-11.2+0.3 


337+10 


-0.04 


06-09-2007 


54350.55 


SKI 


5165-7240 


2200 


-12.6+0.4 


230+10 


0.05 


12-09-2007 


54356.98 


FLW 


6050-7050 


7200 


-12.3+0.5 


275+5 


0.11 


04-10-2007 


54378.49 


SKI 


3265-5310 


1200 








04-10-2007 


54378.53 


SKI 


5165-7240 


2200 


-10.7+0.5 


303+13 


0.06 


25-12-2007 


54460.67 


NOT 


3520-5075 


800 








08-01-2008 


54474.75 


FLW 


4760-6760 


3200 


-6.8+0.2 


324+10 


0.10 


31-01-2008 


54497.78 


FLW 


4780-6780 


3200 


-6.8±0.2 


350+13 


0.08 


02-03-2008 


54528.65 


FLW 


4770-6770 


3200 


-5.8+0.5 


334+25 


0.02 


03-09-2008 


54713.56 


SKI 


5090-7165 


2200 


-2.9+0.1 


370+6 


-0.02 


25-01-2009 


54857.68 


FLW 


4775-6775 


3200 


+ 1.1+0.1 


595+25 


-0.01 


28-09-2009 


55103.57 


SKI 


5020-7100 


2200 


-3.9+0.3 


433+5 


0.0 


14-11-2009 


55150.90 


FLW 


4755-6760 


3200 


-2.9+0.2 


407+5 


0.0 


12-01-2010 


55209.87 


FLW 


4730-6735 


3200 


+0.25+0.15 




0.01 


16-01-2010 


55213.71 


FLW 


4740-6740 


3200 


+0.7+0.2 


525+25 


0.0 


28-03-2010 


55284.85 


NOT 


5785-8285 


1500 


+2.2+0.2 
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Fig. 2. Evolution of the Ha equivalent width (top) and V/R ratio 
(middle) with time. The bottom panel shows the peak separation 
as a function of EW(Hq'). 



a long-term weakening of the disc is suggested by the decrease 
of the equivalent width (Fig.|2]i, which changed from — 12 A in 
September 2007 to — 3 A a year later. Negative values of the 
equivalent width mean that the line displays an emission profile. 
In 2010, EW(Hq') was positive, indicating an absorption profile. 

3.2. Spectral classification 

The blue spectrum of IGR J06074+2205 (Fig.[3] middle) is dom- 
inated by hydrogen and neutral helium absorption lines, clearly 
indicating an early-type star (O or B). The Balmer series lines 
from H/3 up to Hrj at 3835 A are seen in absorption. No or very 
weak He II lines (/14541, /14686) are present, which implies a 
type later than BO. The weak Mg II /14481 indicates a type ear- 
lier than B2. The strong C III -H O II blend allows us to constrain 
better the upper limit as its presence indicates a type earlier than 
B 1 .5. The primary classification criteria in the spectral range 09- 
B1.5, namely the ratio of Si IV at A4686 to Si III at /J4552-68- 
75, cannot be used because the fast rotation blurs the weak Si IV 
line. Nevertheless, Si III 4552-68-75 is clearly present, which if 
the star is on the main sequence, it would imply a spectral type 
earlier than B2; Si III has its maximum strength at type B0.5V 
and is no longer seen at B2V ( Walborn & Fitzpatrick 1990). 

As for the luminosity class, the relative strength of the Si 
III/14552-68-75 complex and the CIII + Oil blends compared to 
that of nearby He I lines favours a main sequence classification. 
Likewise, the low relative intensity of Oil lines also indicates a 
luminosity class V star 

A visual comparison of the I GR J06074-H2205 spec t rum to 
those of the standards shown in I Walborn & Fitzpatrickl d 19901) 
yields that the spectrum of HD 36960 as the most similar one, 
hence we assign a spectral type B0.5V to IGR J06074-I-2205. 
However a slightly later type of Bl V cannot be ruled out. Fig. [3] 
also shows the spectra of two MK standard stars, namely, e Per 
(a BOV star) and 42 Ori (a Bl V star). 



3.3. Reddening and distance 

To estimate the distance, the amount of interstellar extinction 
Ay - Rx E(B - V) to the source has to be determined. The ob- 
served colour of IGR J06074-H2205 is (B - V) ^ 0.58 + 0.03, 
while the expected one for a B0.5V star (B - V)o - -0.26 
(|Johnsonlll9"66l: iGutierrez-Morend [19791: IWegneill l994l) . Thus 
we derive a colour excess of E(B -V) = 0.84 + 0.03. Assuming 
the standard extinction law R - 3.1 and an average absolute 
magnitude for a B0.5V star of My - -3.5 + 0.5 ( Vacca et al] 
1996; Wegner 2006) the distance to IGR J06074+2205 is esti- 
mated to be 4.4+ 1 .0 kpc. Since the photometric magnitudes were 
obtained when the source was displaying a relatively large value 
of the EW(Hq'), this estimate of the distance should be taken as 
a lower limit. In the absence of the disc the colour B - V is ex- 
pected to be smaller (more "blue"), which would imply a lower 
value of E(B - V), hence longer distance. 

The colour excess can also b e estimated from the strength 
of the diffuse interstellar bands dHerbi gl 1 1 975t iHerbig & Lekal 
119911 : iGalazutdinov et alj l2000l) . Using the strongest lines 
(6203 A, 6269 A, 6276-79 A, and 66 13 A) in the red spectrum 
of IGR J06074-H2205, the estimated color excess is E(B- 
V)=0.85+0.12, in excellent agreement with the photometric de- 
rived value, hence giving support to the spectral classification 
suggested. The error is the standard deviation of the measure- 
ments of all employed lines. 



3.4. Rotational velocity 

Be stars are fast rotators. They have, on average, larger observed 
rotational velocities than B stars as a group (Slettebak 1982). The 
determination of the rotational velocity is important because it is 
believed to be a crucial parameter in the formation of the circum- 
stellar disc. A rotational velocity close to the break-up or criti- 
cal velocity (i.e. the velocity at which centrifugal forces balance 
Newtonian gravity) reduces the effective equatorial gravity to the 
extent that weak processes such as gas pressure and/or non-radial 
pulsations may trigger the ejection of photospheric matter with 
sufficient energy and angular momentum to make it spin up into 
a Keplerian disc. At present there is no consensus on how close 
to their critical velocity Be stars rotate, although observations 
suggest that a large fraction of Be stars rotate at 70-80 % of the 
critical value (ISletteba klfigsl iPortedl 1 9961; l^idinlllOOTI) . 

The projected rotational velocity can be estimated by mea- 
suring the full wid th at half maximum (FWHM) of Hel lines 
jSteele et alj|1999l) . The projected rotational velocity was ob- 
tained as the average of the values from three Hel lines at 4026 
A, 4143 A, and 4471 A. We measured the width of these lines 
from the NOT spectrum as it provides the highest resolution 
in our sample. We made three different selections of the con- 
tinuum and fitted Gaussian profiles to these lines. We obtained 
vsin/ = 260 + 20 km s"'. The value quoted is the weighted av- 
erage of the nine measurements. 

Table|3]gives the rotational velocity of some Be/X-ray bina- 
ries obtained from the reference in column 8. Column 6 gives 
the inclination angle of the orbit with respect to the line of 
sight, which allows to estimate the true rotation velocity of the 
Be star The critical velocity depends on the spectral type. For 
Be/X-ray binaries whose spectral type distribution spans a very 
narrow range (09-B2), the critical or break-up velocity is 500- 
600 k m s~', with the lower end corresponding to later spectral 
types (ICranme d2005h . The term shell refers to double-peak lines 
wh ose central depression is lowe r than the stellar continuum (see 
e.g. lHummel & Vranckenll2000h . This type of profile occurs in 
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Fig. 3. Identification of the He and metallic lines used for spectral classification in IGR J06074+2205(iniddle spectrum). A compar- 
ison with two MK standards is provided: the BOV e Per (top) and the BIV 42 Ori (bottom). 



Table 3. Comparison of IGR J06074+2205 with other Be/X-ray binaries. 



X-ray 
source 



Optical 
counterpart 



Spectral 
type 



Disc-loss 
episodes 



Porb 

(days) 



Inclination 
angle (°) 



(km s-') 



Reference 



IGR J06074+2205 




B0.5IV 


yes 






260±20 


this work 


4U 01 15+634 


V635 Cas 


B0.2V 


yes 


24.3 


43 


300+50 


1 


RX J0146.9+6121 


LS I +61 235 


BlIII-V 


no 






200+30 


2 


V 0332+53 


BQCam 


08-9V 


no 


34.2 


<10 


<150 


3 


X-Per 


HD 24534 


09.5III 


yes 


250 


23-30 


215+10 


4,5 


RX J0440.9+4431 


LS V +44 17 


BlIII-V 


yes 






235+15 


6 


lA 0535+262 


HD 245770 


09.7III 


yes 


111 


28-35 


225+10 


7,8 


RX JOS 12.4-3 114 


LS 992 


B0.5III-V 


yes 


81.3 




240+20 


9 


lA 1118-615 


Hen 3-640 


09.5IV 


no 






300 


10 


4U 1145-619 


V801 Cen 


B0.2III 


no 


187 


<45 


250+30, 290 


11,12 


4U 1258-61 


V850 Cen 


B2V 


yes 


132 


90 shell 


<600 


13 


SAX J2103.5+4545 




BOV 


yes 


12.7 




240+20 


14 



[1] Neeueruela et al. (2001), [2] Reig et al. (1997), [3] Negueruela et al. (1999), [4] Lvubimkov et al. (1997) 
[5] Delgado-Marti et al. (2001), [6] Reig et al. (2005), [7] Haigh et al. (2004), [8] Grundstrom et al. (2007) 
[91 Reigetal. (2001) . [10] Janot-Pach eco et al. (1981) . [Ifl Janot-Pacheco et al. (19821 [121 Webster (1974) 
[131Parkes etal(1980), [14] Reig et al. (2004) 



systems with high inclination angles, i.e., the so called edge-on 
systems. 



4. Discussion 

All spectroscopically identified optical companions of Be/X-ray 
binaries in the Milky Way have spectral type earlier than B3. 
To the authors' knowledge there is no one single exception. In 
terms of the mass, this means that there are no Be stars in Be/X 
binaries with masses lower than 8 M©. Thus the classification 



of IGR J06074+2205 as a B8III reported bv lMasetti et aP (l2006l) 
was surprising. 

This narrow range of spectral types in Be stars which are 
part of binary sy stems c ontrast with the wide range found in iso- 
lated systems (iNegueruelal 1 1 99 Si) . The "Be phenomenon" can 
be observed from late O stars to early A stars. Such narrow 
range in masses for Be/X-ray binaries ca n be explained by invok- 
ing a non-conservative binary evolution. Portegi es Zwartl (Il995h 
showed that if during the evolution of the binary the amount 
of angular momentum lost per unit mass through the second 
Lagrangian point increases then the number of late-type Be stars 
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with a neutron star companion reduces. The reason for this is 
that a significant loss of angular momentum implies that a larger 
number of binaries ends tip as mergers. On similar grounds, 
Ivan Bever & VanbeverenI d 19971) compared two models that only 
differed in the amount of angular momentum lost during non- 
conservative Roche Lobe overflow and found that the minumum 
mass of the mass-gainer star corresponding to the model with 
the higher loss of angular momentum was ~ 7 M©. The spec- 
tral type that we determine for the optical counterpart of IGR 
J06074H-2205, namely B0.5V is in good agreement with the 
spectral type distribution of galactic Be/X-ray binaries. 

The Ha line shows V/R variability, that is, the cyclic varia- 
tion of the relative intensity of the blue (V) and red (R) peaks 
in the split profile of the line. The V/R variability is believed 
to be caused by the gradual change of the amount of the emit- 
ting gas approaching the observer and that receding from the ob- 
server due to the precession of a density perturbation in the disc. 
Double-peak symmetric profiles are expected when the high- 
density part is behind or in front of the star, while asymmetric 
profiles are seen when the high-density perturbation is on one 
side of the disc. More precisely, when the high-density part of 
the disc is moving toward the observer we expect to see a blue- 
dominated V > R profile, while when the high-density part is 
receding from the observer red-dominated profiles V < R are 
expected ('Telting et al.' 1994). For systems with high inclination 
angles, the two symmetric cases can be readily distinguished 
since the central depression between the two symmetric peaks 
is much more pronounced (reaching or going beyond the contin- 
uum), adopting a shell profile, when the perturbation is behind 
the star If the density perturbation revolves around the star in the 
same direction as the material in t he disc (prograde pr ecession) 
then the V/R sequence would be dTelting et al]lI994 : V - R 
(perturbation behind the star) — > V > R — > V - R (perturba- 
tion in front of the star, shell profile) — > V < R. 

Figure [T] shows the Ha line profile from most of the spec- 
tra that we obtained. The first recorded spectrum displays a red- 
dominated peak, V < R. The relative intensity of the two peaks 
appeared reversed, i.e. V > R, about a year later If we inter- 
preted the January 2009 line as an emission line affected by a 
strong shell component (self-absorption from the disc), then the 
V > R would be followed by a shell phase and we would con- 
clude that the motion is prograde. 

The available data do not allow us to constrain the duration 
of the V/R cycle because a major structural change occurred 
which brought the V/R cycle to an end. This major event was 
the loss of the equatorial disc. Nevertheless, we can constrain the 
quasiperiod of the V/R cycle from the two complete phases cov- 
ered by the data, namely, the V > R and shell phases. The shell 
phase would have lasted for about 14-16 months (July- August 
2008 to November -December 2009). Due to the observational 
gaps the V > R phase is less constrained. The maximum dura- 
tion for this phase is approximately 20 months (December 2006 
to August 2008). Since we do not expect the transition between 
phases to be abrupt, it reasonable to limit the V > R phase to 15- 
17 months (February 2007 to June 2008). Assuming then that 
each V/R phase lasts for about 15 months the V/R quasiperiod is 
estimated to be ~ 5 years. 

The long-term decline of the Ha equivalent width (Fig. [2]i 
undoubtedly implies a progressive weakening of the disc over 
a period of ~3 years. Moreover, the latest spectrum in March 
2010 displays an absorption profile. The reversion from emis- 
sion to absorption is generally believed to be due to the dis- 
appearance of the disc. Nevertheless the fact that the line still 
shows some residual emission at the wings and that, although 



positive, the Ha equivalent width, +2.2, is still somehow smaller 
than expected for a non-emitting B0.5 V star, which according to 
Jaschek & Jaschek ( 1987) should be ~ +3.5-4 A may imply that 
the disc-loss phase continues at the time of writing. 

Irrespectively of whether the disc vanished completely or 
not, the data show one important result, namely, the V/R vari- 
ability is absent in small discs. In other words, only in well 
developed discs can the density perturbation survive. We have 
searched the literature for spectral line variability immediately 
before and after disc-loss episodes in other Be/X-ray binaries 
(see column 4 in Table [3]). We observe that shortly before the 
diappearance of the disc and during the first instances of the 
formation of a new disc the Ha line profile is symmetric and 
double peaked. None of the Be/X-ray binaries that have gone 
through disc-loss phases and for which there is a good op- 
tical follow-up coverage, namely, X Per in 1990 (Clark etal] 
2001j), 4U 0115+63 in 1997 lNegueruelaetalJl200ll) and in 
2002 (iReig et al.' '2007b lA 0535+262 in 1999 (iHaigh et al.1 
1999; Grundstromet aLl l2007h and RX J0440.9+4431 in 2001 
(Reig et al. 2004) exhibited asymmetric profiles during the ini- 
tial stages of disc formation. That is, the effects of the density 
perturbation do not show up until the disc is fully developed. 
The same result holds for the final stages of disc loss. Just be- 
fore the detection of the absorption profile the V/R ratio is ~ 1 . 

Using the equation given in Sect. 13.11 the rotational veloci- 
ties given in Table [3] and the peak separation of the first avail- 
able asymmetric profile after the disc loss Apeak = 220 km s ' 
(X-Per, Clark etal."200l'). Apeak = 305 km s"' (4U 0115+63, 
Negueruelaet al. 2001) we find that the disc must reach a ra- 
dius of ~ 4Rt for the onset of the V/R variability. Virtually the 
same value is derived from Fig. 5 in'Grundstrom et al.' ("2007^ for 
lA 0535+262. For RX J0440.9+4431 (Reig et al. 2004) only a 
lower limit, > 2/?,, can be given. 

The data before the disc-loss phase is even more scarce. 
Nevertheless, the available data seems to indicate that, before 
the disc loss, the density wave may survive in smaller discs 
(~ 1.5 - 2Rt). In 4U 01 15+63, the last available spectrum prior 
to the disc-loss event showing substantial V/R variability corre- 
sponds to a disc radius of 1.6 In IGR J06074+2205, the peak 
separation measured in the spectrum taken on 1 January 2008 
corresponds to radius of 2.6 R^- 

Note that the peak separation does not necessarily indicate 
the geometrical size of the disk, but rather the radial distance at 
which Ha emission becomes optically thick. How well this cor- 
relates with the disc size depends on its density distribution. It is 
generally assumed that the density falls off as a power law with 

increasing distance from the star, p(r) — po (^^^ . The fact that 
the power-law index n est imated obseryationally lies in a very 
narrow range 3 < n < 3.5 (Water s et all 1988t Hones et al.ll2008f) 
and that the relevant radius in this order of magnitude calcula- 
tion is that of the Ha emitting region justifies the validity of the 
comparison. 

Given the large observational gaps in the present observa- 
tions it is clear that further long-term spectroscopic monitoring 
of these and other systems are required to confirm these results 
and set better constraints on the variability timescales and be- 
haviour of the disc perturbations. 

5. Conclusion 

We have performed optical photometric and spectroscopic ob- 
servations of the optical counterpart to IGR J06074+2205. From 
the photometric magnitudes and colours and the strength of var- 



Reig et al.: The optical counterpart to IGR J06074+2205 



7 



ious diffuse interstellar bands we have estimated the distance to 
be ~ 4.5 kpc. From the ratios of various metallic hnes we have 
derived a spectral type B0.5V and from the width of three He I 
lines estimated the rotation velocity of the underlying B star in 
260 km s-'. 

The long-term optical variability of this system is charac- 
terised by global changes in the structure of the equatorial disc 
around the Be star companion. These global changes manifest as 
asymmetric profiles of the Ha line and a significant decay of its 
intensity. The asymmetry consists of V/R variability with char- 
acteristic timescales of the order of few months. The Ha line 
has been seen to change from emission to absorption, indicating 
the loss of the equatorial disc. The V/R cycle stopped when the 
equatorial disc became too small. A small disc cannot support a 
density perturbation. During the final stages of the disc evolution 
V/R asymmetries can be seen up to a disc radius of ~2 The 
density wave fades away before the complete dissipation of the 
disc because the disc became too tenuous to support the density 
perturbation. A comparison with other systems reveal that, after 
the loss of the disc, a density perturbation does not develop until 
the disc radius is of the order of 4 
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